CRE components is overestimated and the other underestimated, both cannot be driven towards 
observed vaiues by adjustments to cloud condensate heterogeneity and overlap alone. 
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Abstract 


The radiative impacts of introducing horizontal heterogeneity oftay^cloud condensate, and 
vertical overlap of condensate and cloud fraction aic ovuvned with tMutipn a new radiation 
package operating in the GEOS-5 Atmospheric^retieraMTtiCi’lation l^a'del. The impacts are 
examined in terms of diagnostic top-of-tl^:»tniwp!ie^shottwB c^SW) and longwave (LW) 
cloud radiative effect (CRE) calculations for a tftgeal assumprtions and parameter 


cloud radiative effect (CRE) calcuMtfhs for a assumprtions and parameter 

specifications about the overlap. 3Tt^/,im^estigation wBonducted for two distinct cloud 
schemes, the s^^MLGEOS-5 distribution, and another which has 


been recently - Wed experimentally frarats cnhaMm’ cloud microphysical capabilities; both are 


coupled d^ud generator allowing arbitrary cloud overlap specification. We find that 
cloud overlapnrawfoye impacts are sLujficantly stronger for the operational cloud scheme for 
which a change oMjoiid fraction »srlap from maximum-random to generalized results to 
global changes of SW JUd of ~4 Wm' 2 , and zonal changes of up to ~10 Wm' 2 . This 

is because of fewer occurT'.'jJKs compared to the other scheme of large layer cloud fractions 
and of multi-layer situations with large numbers of atmospheric being simultaneously cloudy, 
conditions that make overlap details more important. The impact on CRE of the details of 
condensate distribution overlap is much weaker. Once generalized overlap is adopted, both 
cloud schemes are only modestly sensitive to the exact values of the overlap parameters. We 
also find that if one of the CRE components is overestimated and the other underestimated, 
both cannot be driven towards observed values by adjustments to cloud condensate 
heterogeneity and overlap alone. 
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2 1 Introduction 

3 With new computationally efficient approaches to treat cloud-radiation interactions now 

4 available, there are fewer reasons to retain the simplistic cloud descriptions that have persisted 

5 in General Circulation Models (GCMs) for many years. Clouds do no longer have to be 

6 treated by the radiation schemes of these models as homogeneous slabs within large areas 

7 O(10 4 km 2 ), with fractional coverages and optical depths that have been greatly adjusted to 

8 compensate for known biases arising from their nonlinear HEcftctiou wan radiation. While 

9 capturing the radiative effects of full-blown 3D cloud hwrogeneitwRwtill be elusive, the 

10 representation of in-cloud horizontal heterogeneity ofTCloud^ondeny® f md tVo-point 

11 statistics of vertical correlations of condensate and cloud i«u within a on|rdu«' ensional 

12 radiative transfer framework is now feasible. As a matter of fact, t?$ current work is one more 

13 study that amply demonstrates the viability of suefrjia midkrtakmg. <£ 

14 The main development that makes more ootoplcx de&rijftions possible is the 

15 introduction of methods that perform rajiiljW nan-stm in thi vluady portions of GCM grid 

16 cells in a stochastic manner (Pincus c( 21)03). I h^mnre.ct'mplcx cloud descriptions come 

17 from cloud generator^ priMhu^u h^riz^rtM^nd vertical cloud variability according to rules 

18 that are relatively -easy touin»lenNwU. Th^blbud/ields from the generators can then be 

19 coupled wit]* Stochastically o perarfl«ud iatt ve 1 ransf er schemes that only “see” atmospheric 

20 subcolqjjra^^^ne cloud fraction is limy and condensate is horizontally invariable whenever 

21 a layer is cloudy Wjth the radiative ttdnsfer simplified, the sensitivity of the radiation budget 

22 to a variety of speci.’k uions tharaansform a gridcolumn’s cloud profile to a cloud field 

23 consisting of several sraRjdjfms can be easily examined. What we should ultimately 

24 investigate is whether thejsifects of cloud complexity on the transfer of solar and thermal 

25 infrared radiation matter for the GCM’s climate. Such a study of the full impacts of 

26 interactions and feedbacks of the altered radiation fields with the multitude of the GCM’s 

27 dynamical and physical processes is left for the future. Here, we simply focus on diagnosing 

28 the possible range of radiative impacts of enhanced cloud complexity, an approach akin to that 

29 of Shonk and Hogan (20 1 0). 

30 In the following we will present the tools, assumptions, and experimental setup that allow us 

3 1 to examine the degree to which cloud complexity changes the cloud radiative impact (sections 
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2, 3, and 4). The availability of two cloud schemes in the GCM at hand and our analysis 
approach provides the opportunity to investigate whether the same assumptions about cloud 
complexity imposed on different initial cloud fields can yield notably distinct radiative 
impacts (section 5) and what causes the contrasting behaviour (section 6). 


2 Implementation of RRTMG into GEOS-5 


The effects of cloud overlap (fraction and condensate) on the radiative fluxes can be captured 
best with radiation codes equipped with as much flexibility in tbe representation of 

such overlap. This (along with improved representation oagaseous afcxyknon) was one of the 
primary motivations for the implementation into tl$»£]EQ^-5 ^Boitpheric _ General 
Circulation Model (AGCM, Rienecker at al. 2010; Mo^ b al. 2012) 'V Ifr^fLRTMG 
radiation package (Clough et al 2005), a faster version Or thcN^.R|^/f codes (Mlawer et al 
1997; Iacono et al. 2008) designed specifically for lj 
and thermal infrared component. Both comp 
Independent Column Approximation (McJT’A 1 


MclCA has been implemented succesfthfT into EC] 

/ ' * Jr, 

(Morcrette et al. 2008) and several ^Ihet^rge scale 
when the radiation on .V j 

integration ovosSeach column is'^pj^fed by s 



;ale m mids^d ^pjysisting of a solar 


.^tfcnression of this process 1 can be -iVritten as follows: 


run in sA' 'called Monte Carlo 
al. 2003). RRTMG with 
ntegrated Forecasting System 
Within the MclCA framework, 
ber of atmospheric (sub)columns, full spectral 
ic (Monte Carlo) integration. A simplified 



n \ 




n«l W ' w-U-1 k~ 1 

The uppercase symbols represent broadband fluxes, while the lowercase letters 

represent pseudo-monochromatic fluxes per the correlated-k paradigm (Lacis and Oinas, 
1991). F represents a broadband flux (solar or thermal infrared; upward or downward) at any 
vertical level within the AGCM gridcolumn, F„ is a similar broadband flux for one of the N 
subcolumns generated by RRTMG’s cloud generator (Raisanen et al. 2004, see below) within 
the gridcolumn, and/„ >(t is the pseudo-monochromatic flux for subcolumn n and spectral point 
k. What the above equation essentialy conveys is that a broadband flux which is normally 
obtained by taking the average over N subcolumns of the sum of K spectral calculations for 
each subcolumn, is approximated by the stun of K spectral calculations where each spectral 
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1 point k is paired randomly with one of the N subcolumns, n*. Note that when using eq. (1) the 

2 computational cost of the calculation over all subcolumns is the same as that of a full spectral 

3 integration of a single (sub)column. The performance of this approximation in large scale 

4 models has been tested extensively (e.g., Barker et al. 2008). The main issue of concern is 

5 whether the conditional random noise, decreasing as the inverse square root of the number of 

6 applications of eq. (1), has any detrimental impact on the simulations. The prior studies and 

7 our test with GEOS-5 have shown that the MclCA noise for sufficiently long runs (at least a 

8 month) is of similar length and nature as the internal variability of the model. 

9 An extensive description of the brand of generator used ir. otilSmp 1 gestation of RRTMG is 

10 provided by Raisanen et al. (2004). The cloud generator Educes sq^Wturvj’ i s that have either 

1 1 clear or completely overcast cloud layers. Whether the cloud |*uensate or <t (vrticiilar layer 

12 is different from one subcolumn to the next depends onjplw,v3>;itnptions atwjflPmrizontal 

13 cloud heterogeneity: either homogeneous or heterogeneous coru^Rtac distributions can be 

14 specified within the cloud generator. The horizontal loclujon of clouS, 1 . uPa. particular layer 

15 (i.e., subcolumn assignment) and specific valueAf coni^muic^fordieWrogeneous condensate 

16 distributions) depend on cloud presencjfrta omer fryers scoping to the overlap rules 

17 implemented. By design, in the limiltejfan infinite fimhj^fcf subcolumns layer horizontal 

18 averages reproduce the vcfiiuql nrow1e^Bck>ud fraction ami condensate provided as input to 

19 the generator by £jf»> SG^^Nfri^^DecificWt&iptions of rules and assumptions about cloud 

20 fraction and ^undunsate distribution Burlaps asrilfiplemented in the GEOS-5 cloud generator 

2 1 are provided in fee section that folio* « 

22 if 

23 3 Cloud overlap and variafckity representation 

24 The cloud fraction overla p, options for the cloud generator that comes with the RRTMG 

25 package include the standard assumptions that have been used extensively in the past, i.e., 

26 maximum, random, and (the most popular) maximum-random overlap (Gcleyn and 

27 Hollingsworth 1979; Tian and Curry 1989) where contiguous cloudy layers overlap 

28 maximally and randomly otherwise. Raisanen et al. (2004) provides a mathematical 

29 description of the practical implementation of these overlap assumptions in a cloud generator 

30 algorithm. In this work, from the above simplified overlap descriptions, we only test the 

3 1 maximum-random overlap option. 
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Starting with the work of Hogan and Illingworth (2000), numerous studies (e.g., Mace and 
Benson-Troth 2002; Oreopoulos and Khairoutdinov 2003; Naud et al. 2008) have shown that 
the above simple overlap assumptions are inconsistent with cloud fields from observations 
and cloud resolving models, and that the concept of “generalized” cloud fraction overlap 
represents observed overlap more realistically. In the generalized overlap paradigm, the 
combined cloud fraction of two cloudy layers at heights z\ and Z 2 with separation distance A z 
= Z 2 -Z 1 can be approximated as a weighted average of combined cloud fractions from 
maximum and random overlap, C max (Az) and C ron (Az), respectively according to: 


C(Az) = a(Az)C mijr (Az) + (l - a(Az))C ra „(Az) 

with 

C mrat (Az) = max(C(z 1 ),C(z 2 )) 

C„„(Az)=l-(l-C(z l ))(l-C(z 2 )) 

The weighting parameter a(Az), is a mea^jr dt'djjc pflxirni^ u< overlap to maximum (exact 
when a(Az) = l) or random (exact whs" ’rMz) ~G); KCgativi. values suggest some degree of 
minimum overlap (a combined cyHnWction greai«®Km that of random overlap). A 
commonly used att^liBwuiOfvBWso aoB| ! here, is that a(Az) depends only on the 
separation disttmce A z and not oi^K^ptcific^Spcs of z\ and Z 2 , i.e., cloud fraction overlaps 
the exact <i:iievv:iy at different height* of the atmosphere as long as A z is the same. With this 
assumption, ii^ir. \hown (Hogan md Illingworth, 2000) that a(Az) can be fit reasonably 
well by an inverse ^“Jfev^ential tuncflpn: 



a( Az) = exp - 


Az 


( 4 ) 


22 where L a is the “decorrelation length” for cloud fraction overlap. Such a fit obviously does 

23 not allow for negative values a(Az) which are occassionally observed (e.g., Oreopoulos and 

24 Norris 201 1). Because the fit provided by eq. (4) is usually used in conjunction with eq. (2), 

25 generalized overlap has also been termed “exponential-random” overlap (Hogan and 

26 Illingworth 2000). 
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The manner in which cloud water contents align in the vertical may also important for 
processes like radiation (or precipitation). For example, the domain-averaged fluxes differ 
between a case where all high or low condensate values are aligned to create pockets of 
vertically integrated high or low liquid water path (WP), and a case where a more random 
alignment homogenizes the WP horizontal distribution (e.g., see Norris et al., 2008). The 
nature of condensate alignment can be expressed in terms of rank correlations of water 
content as a function of separation distance A z = zj-z\ (e.g., see Pincus et al., 2005 and 
Oreopoulos and Norris 2011). For two layers at heights z\ and Z 2 the water contents at both 
heights can be ranked separately for the overlapping portioi^fcA &ks ufo olumns of the two 
cloud layers. A linear correlation coefficient r(Az) can the* rm the ranks Ri(z\) 

and Rfa) according to: 


r(Az) = • 


u 

x=l 




The rank correlation coefficient expressc - ™ likelih&d wa^r contents of the same relative 
strength within their respective layer .ira Aligned in uiKvcrtic il, with r(z\,Z 2 )=\ corresponding 
to perfect alignment to completely random alignment. 




It was suggesidfl (e.g.. RSisanen o'lO 2004) ihiflflilk; rank correlation coefficient can also be fit 

V ? A Jr 

by an invciiC- exponential (which .igliin will not capture negative values) under similar 

assumptions al- ijfjlihe case of the eh fid fraction overlap parameter, i.e., that it is only a 

m 

function of A z and M! i'itself 




r( Az) = expj - - 


( 6 ) 


where L r is the rank correlation decorrelation length. Large values of L r indicate condensate 
values that are highly correlated in terms of relative strength, while small values suggest 
condensate values whose relative strength is weakly correlated between layers. 

The practical implementation of generalized cloud fraction overlap and condensate overlap 
using inverse exponential fits is described by Raisanen et al. (2004). The cloud generator that 
came with RRTMG could handle generalized cloud fraction overlap, but did not allow for 
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overlap of condensate distributions; we added that feature following Rdisanen et al. (2004). 
To create the subcolumns that describe the cloud fields within the GCM gridcolumns, two 
additional pieces of information, besides the profiles of cloud fraction C and mean condensate 
(liquid and ice) are needed, namely to specify the decorrelation lengths L a and L r and the 

magnitude of the horizontal variability of the condensate distributions. We defer discussion of 
decorrelation lengths for the next section, and describe variability here. 


To create condensate distributions for cloudy layers we assume that beta distributions describe 
the horizontal variations of normalized condensate x=wlw max \ 


v (x)= RP±3l x r - i 

Ppi T(p)T(q) 


(l-x) 


i - 1 




( 7 ) 


where T is the gamma function and the maximum value of c.'t^igj^ate w max is set .r< five times 
the assumed variance of the distribution. The shape parameter;- ip ...of the beta distribution 

are calculated from the method of moments (Wilks 


x 2 (l - jc) 


where x = w/n 



(8a) 

(8b) 


The standard tfey uition a w of the diiifthlrtion was set as follows, loosely based on Oreopoulos 
and Barker (13W'( and our own anal y J of hydrometeor variability in the CloudSat (Stephens 
et al. 2002) data: 

a w =0.5'! vvtii’i. ?/ > 0.99 

o w = wl*j2 .'.hen 0.9 sCs 0.99 (9) 

ol = w when C < 0.9 


The choice of the beta distribution supported by observations (Oreopoulos and Davies 1998 
and Lee et al. 2010), but other skewed distributions that have also been observed from 
airborne and satelllite measurements, such as gamma and lognormal would have been an 
equally acceptable alternate choice. Eqs. (7) and (9) apply to both liquid and ice condensate, 
and in layers where the two phases coexist their ratio is assumed to remain constant across all 
subcolumns. Since no distinction is made between liquid and ice cloud fraction, the 
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normalized standard deviation a w /w is de facto the same for liquid and ice condensate 
distributions. The beta distribution of normalized condensate x is converted to an actual 
condensate distribution and then to a cloud optical depth distribution using the AGCM- 
provided effective particle size which is different for each phase, but assumed horizontally 
homogeneous. 


Since the specification of the amount of condensate variability via o w does not come explicitly 
from the host AGCM or derived from rigorous physical principles, and variability is used only 
to gauge diagnostically the sensitivity of the cloud radiative effect, we argue that it is not 
essential to fully justify its exact specification. Di ffereB^SJffecJ of .Variability will have 
quantitatively different impacts on the cloud radiative Ijfect, butane! tjWhtative impact is 
nevertheless entirely predictable: larger inhomogeneity res^iji'ili Smaller sMlvfave (SW) and 
longwave (LW) cloud radiative effect and vice-versa. . , , 

4 Description of AGCM setup and expeVtrrscntS Jr 


4.1 Specification of overlap parameter decorrelation lengths 

As explained earlier, for. the ■■rrients with generalized cloud fraction overlap and 

heterogeneous ciX^trllsafe' di.Mi ihhkons, tnfc decorrelation lengths L a and L r need to be 

specified. t* ■rnplest option is ttfjsclect valucJWiat are universal (constant) in space in time. 

reen used in prioiflumk (Raisanen et al. 2004; Morcrette et al. 2008) are 
andT®L km. Such a.smplification is not justifiable in principle on physical 

k w 

grounds given the wuJ, range -.of cloud regimes. Still, whether a more sophisticated 
specification of decorrelaiiniu .lengths is needed in practice requires further investigation. The 
availability of cloud partjcle/hydrometeor reflectivity and backscatter data from the Cloud 
Profiling Radar (CPR) of the CloudSat mission and the CALIOP lidar of the CALIPSO 
mission (Winker et al. 2010), potentially allows a more detailed examination of 
spatiotemporal variation of cloud overlap decorrelation lengths. 


Values ihai 


L a = 2 km 


We performed such a cloud overlap analysis using CloudSat products for two months, January 
and July 2009. For cloud fraction overlap we used the 2B-GEOPROF-LIDAR product which 
provides a cloud mask from combining the different hydrometeor detection capabilities of 
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1 CPR and CALIOP (CPR is more capable at detecting layers with large concentrations of 

2 hydrometeors while CALIOP can better detect unobscured thin clouds). For condensate 

3 distribution overlap we use CloudSat’s 2B-GEOPROF product which provides reflectivities 

4 for footprints (~ 1.7 km) that have been identified to contain hydrometeors at various vertical 

5 locations (separated by ~500 m). Our rank correlations according to eq. (5) therefore actually 

6 come from reflectivities and not cloud condensates which are also available from CloudSat 

7 (e.g. product 2B-CWC-RO or 2B-CWC-RVOD), but are considered less reliable for the liquid 

8 phase due to drizzle and mixed/supercooled clouds often assigned erroneously to the ice phase 

9 
10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

22 (~170 km), similar die spatial resqijliion of the AGCM experiments described below. There 

23 is a clear zonal strucSb^&t both -ninths with tropical latitudes exhibiting larger decorrelation 

24 lengths (more maximum ovrfpap and greater vertical alignment of reflectivities of similar 

25 relative strength), consisted' with overlap contrasts between convective and stratiform regimes 

26 (Barker 2008a, b; Oreopoulos and Norris 2011). L r values seem to be generally about half 

27 those of L a , in broad agreement with previous findings (Raisanen et al. 2004; Pincus et al. 

28 2005; Oreopoulos and Norris 201 1). Seasonal shifts of the peak values of decorrelation length 

29 appear to reflect the movement of the Intertropical Convergence Zone (ITCZ). 


(Lee et al. 2010). Since reflectivies are proportional to the-jtj/o »)f hy^ometeor particles, 
under the assumption of constant particle number, the amtunt of coadnmte is monotonically 
related to particle size and eq. (5) can be applied to reflec®kie^^ well. A £fcveat of the 2B- 
GEOPROF reflectivities on the other hand is that they do nut Jesuit only frou^ftwittions of 
the radar beam with suspended particles, but also precipi tatiofl^Wi tides. While the above 
make CloudSat-derived decorrelation lengths apprc Aini*^, it should^^feic/in mind that the 


goal is not to obtain a perfect map of their geog 


plausible general picture of their spatial sd 
globally constant decorrelation lengths^rjfSTb u d radia^ ^g effg vt studies. 

Fig. 1 shows the zona 
square fits (Pres* 
profiles wit 


m, but to liave at our disposal a 
at can be contrasted with 


' ' panel) and L r (bottom panel) derived via least- 

^ontlft^K^f^ifd CloudSat'CALIPSO a(Az) and r(Az) 
Januai ',.rcind July (solid lines), with the limitations 


expIain^fUeurfyer >Jhe data segment Icnlph used in the above calculation is 100 CPR profiles 



\ , 

3 C latitude zones, 
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1 Our objective for AGCM parameterization purposes is to capture the observed decorrelation 

2 length zonal structure shown in Fig. 1. For that purpose, we apply a Gaussian fit (black 

3 dashed curves) of the form 

4 L = m, + m 2 exp[-(0 -m-^/m*] (10) 

5 to the January (black) curves. In eq. (10), 6 is the latitude in degrees and mi, m 2 , m3 and W4 are 

6 parameter fits. All, except m3, are held constant, and their values yielding decorrelation length 

7 in km are provided in Table 1. Parameter m3, regulating the lati tude at which eq. (10) peaks, 

8 reflects the zonal movement seen in the CloudSat data, an4^ anowtM hy vary as a function of 

9 the day of the year according to: 

10 m 3 = -4m i0 (jday -272)/365 when jday >181 '\ 

11 m 3 “ 4m 30 (y£&zy -91)/365 whenjday sl81 ' -. (lib) 

12 where jday is the julian day. We set m3 ( o®.0 (clffllU faction overlap) and m3,o=8.5 

13 (condensate/reflectivity overlap). Our aprwwclw! essence of assigning the initial 

14 Gaussian fit of the monthly-averaged &r»ary observ Bftfis ^ January 1, and then finding the 

15 zonally-averaged decorrelations for days of tn^HJr by applying eqs. (10) and (1 1). 

16 This is the way curves 1 (for July 1) were derived. Note that the 

17 January fits defcribe the zonal a^tjmttion oHb^n decorrelation lengths more realistically 

18 than the July Curves which are no I fit; to the data, but outcomes of the parameterization 

19 expressed by eB 1 1 0) and (11); tlilc parameterized northward shift of the January curves 

20 intended to capture Juiy^verlap gewally leads to underestimates. Again, for the purposes of 

21 this study, where the gc-al W:,^!>&rnine the sensitivity of the cloud radiative effect to a range 

22 of decorrelation length specifications and the differences arising when the exact same overlap 

23 assumptions are applied fo two different cloud schemes, we consider imperfect matching to 

24 observed overlap (itself coming with its own limitations) acceptable. 

25 4.2 Description of AGCM experiments with diagnostic radiation 

26 To examine the changes in the radiative impact of clouds when different assumptions are 

27 invoked about (a) the horizontal heterogeneity of their condensate; (b) the way their 

28 condensate distributions overlap; and (c) the way their cloud fractions overlap, relatively short 
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1 (~1 year) simulations with the GEOS-5 AGCM are conducted with RRTMG producing 

2 “diagnostic” only fluxes. Had we wanted to examine the full impact of our cloud changes on 

3 the model climate much longer simulations of at least a decade with interactive RRTMG 

4 would have been necessary. By diagnostic RRTMG radiation fields we mean that the heating 

5 and cooling rates produced by the RRTMG calculations are not supplied back to the AGCM 

6 to affect dynamical and physical processes. Instead, the model run is driven by the radiation 

7 fields of the original (operational) radiation package (Chou and Suarez 1999; Chou et al. 

8 2001) which treats clouds according to its default configuration, as usual. The McICA 

9 configuration of RRTMG simply runs side-by-side with dagffiffi’.iha.l radiation package and 

10 operates on the cloud fields produced by the standard mcgel. but asMOTBfcrmed by the cloud 

1 1 generator in accordance with our heterogeneity and overlap Hy>unpraons. \ 

N ,;v \, 

12 Our suite of experiments is summarized in Table 2. All exp cna^mt were run vu)i jPe GEOS- 

13 5 AGCM Fortuna 2.5 at 2x2.5° resolution with 72 vertical lev^, feud differ only in their 

14 assumptions about cloud fields. While all exppsjmcntii'^liare the stme profiles of cloud 

15 fraction and mean condensate, other assumption^ itboutflife nature of tie clouds are different 

m 

16 from experiment to experiment. Clouds rffiuiwd to t*J;:>rizontally homogeneous or 

17 heterogeneous and their cloud fractiovtit^aii overlap awbuilfip to either the maximum-random 

1 8 or generalized nverlan napalms. Wnwn wupds are heterogeneous and overlap according to 

19 the maximum-raurfo^o^rhttr in jn t i on/tHghid ensate decorrelation length still needs to be 

20 supplied. Al) >: emulations correspond 1 ^) B-mrKffTruns from which the last 12 months are 

21 considered h’Jalu.i lysis; prescribed t urface^ temperatures for the period May 1993 to May 

22 1994 are used. \ 

23 Two sets of experinKntvyvcie cv.jjuucted. One where the standard (control) cloud scheme 

24 (Molod et al. 2012) oper^Sfliumone with McRAS-AC (Sud et al. 2012, in preparation; Sud 

25 and Lee 2007). The two £Kiud schemes share the same convective scheme (RAS), but with 

26 different assumptions about the onset of convection, and ambient air entrainment (quadratic in 

27 McRAS versus linear in standard RAS) and are fundamentally different in their stratiform 

28 cloud parameterizations and microphysics descriptions. The control cloud scheme has pre- 

29 specified liquid and ice particle sizes, while McRAS-AC has active cloud microphysics where 

30 condensate amounts, particle sizes, and precipitation depend on the aerosol loading. For our 

3 1 experiments we chose to provide McRAS-AC with a present day climatology of aerosol mass 

32 concentrations produced by the GOCART (Chin et al. 2000) chemical transport model. Note 
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that for both sets of experiments, while the aerosols are radiative active in the operational 
radiation package that provides interactive radiation fields, they are not considered in 
RRTMG which produces the diagnostic radiation fields used to assess overlap radiative 
impacts on CRE. 


For each of the experiments we generate the monthly, seasonal and annual geographical 
distribution of the longwave (LW) and shortwave (SW) cloud radiative effect (CRE) at the top 
of the atmosphere (TOA). The CRE is defined as: 


CRF = F ch - F cld 

^ lxlj LW.SW 1 LW.SW 1 I.W.SW 


(12a) 


which can also be written as 


CRE 


l.W.SW 


dr r* ovc \ 

lotWLWSW r LWj5WJ 




liiw 


where F is the outgoing flux (LW or SW) at the TOAj't,> designate \w^(£loudless) skies, 
cld a mixture of clear and cloudy skies, and oWpvercSftjgs (100% Jroud fraction); C,„, is 
the total vertically projected cloud fraction CW ft witty s comes from eq. (12a); 

nevertheless, eq. (12b) which applies wS^nhc clouAsky tfiix can be written as the linear 
combination of clear and overcast can be ulfed ft>f interpreting the CRE, since a 

gridcolumn’s C, ot qguut uniquciy.--de fined tin; rather depends on the cloud fraction overlap 
assumption (fo 4 ,' die same clouoWdKtt protflfflhe closer the overlap to random, the larger 
Cj). F o^Lt*; co mplete intercompy ivifrn of CRE among all experiments we use globally- 
averaged valthsagor select experiments we also compare zonal (latitudinal) averages and 
geographical distntwutes. Although pot important for understanding the sensitivity of CRE to 
cloud heterogeneity ahd also include in our comparison TOA CRE from the 

CERES EBAF v. 2.6 data jeffLoeb et al. 2009) for the period March 2000 to June 201 1. 


5 Analysis of Cloud Radiative Effect dependencies 


24 5.1 Global changes in CRE 

25 We first focus on the sensitivity of globally-averaged CRE to different assumptions about 

26 how cloud fields can be generated from profiles of cloud fraction and mean condensate. Fig. 2 

1 Minimum overlap of various degrees produces even larger C,„„ but there is no such overlap in our experiments. 
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and Fig. 3 chart this sensitivity for the control (CTL) and McRAS-AC cloud schemes, 
respectively. The center box contains AGCM results for the “default” (reference) 
configuration, namely homogeneous condensate distributions and maximum-random cloud 
fraction overlap (Exp. 1, see Table 2). Blue numbers depict CRElw and red CREsw values. 
This box also contains the observed global CREs according to the CERES EBAF (Loeb et al., 
2009) product. The numbers in parentheses show the difference between the AGCM with 
default cloud configuration and observed CRE (reference model-observations). Because 
CRElw is a positive quantity, model underestimates yield negative differences, while the 
opposite is true for CREsw which is a negative quantity. Th<? jjy %X 9 S dhow how much the 
CRE differs from the reference model results in the cenftr box wiitft^^knmptions about the 
nature of the cloud fields change, i.e., reference (-run V •; .feumogpii o u.v) mi pus other 

experiments. For example, when keeping the cloud fractioj) fl^tcjap maximKjn-rniidom, but 
allowing the clouds to be inhomogeneous according' to eqs ,(7)-(9) (leftmost box, 


corresponding to Exp. 2), CRElw decreases by 2.3 W 
absolute value (i.e., a smaller negative value Vjpv 5.6 
homogeneous but changing tire cloud fraetfr 
L a = 2 km, Exp. 3), CRElw increase Py 4,3 W 


absolute terms by 4.3 
adopted for re 
LW, negative- 
numbers 

effect) appear positive number, i 
interpretation chan); 




rahpp&lso decreases in 
(-5,6). £/hen keeping clouds 
d (with globally constant 
and CREsw also increases in 


Therefore, because of the conventions we have 
the the-, CRE arising from eq. (12) (positive for 

Ei.w (srrpjpger LW radiative effect) appear as negative 
5, while increases in CREsw (stronger SW radiative 
Tien the sign of the differences is reversed, the 
i.e., positive CRElw differences signify weaker LW 


radiative effect, while 


■sw differences also signify weaker SW radiative effect. 


Having clarified the sign ^cuiventions of our CRE differences, we now proceed to the physical 
interpretation of the results. We start with Fig. 2 which refers to the CTL cloud scheme. 
Introducing heterogeneity (inhomogeneity) in the condensate distributions following eqs. (7)- 
(9) reduces the strength of CRE (box 2, corresponding to Exp. 2 on the left). This is because 
for the same mean condensate, heterogeneous clouds reflect less solar radiation (e.g., Cahalan 
et al. 1994) and emit less (transmit more) LW radiation (Barker and Wielicki, 1997). For this 
particular case therefore changes in CRE can be attributed to changes in E^ sw in eq. (12b): 
the SW outgoing flux for overcast conditions goes down, while the LW outgoing radiation 
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goes up; in both cases the contrast with the clear-sky flux is reduced. The change in CREsw is 
more than double that on CREiw since the nonlinearity of the LW emittance curve is restricted 
to a much narrower range of cloud condensates (or, strictly speaking, optical depth) than the 
nonlinearity of the SW albedo curve. In other words, changes in the details of an optical depth 
distribution begin to matter less (because of saturation in emittance) at lower values of mean 
cloud optical depth. When cloud distributions remain homogeneous, on the other hand, but 
cloud fraction overlap changes (transition from box 1 to box 3), it is C to t in eq. (12b) that is 
mainly affected (it appears from our results that the change in the distribution of cloud tops 
exposed to space, which matters for the LW, is a lesser comfit r Mtso t >■ CRE s w and CRElw 
become stronger by the same magnitude (4.3 Wm" ), indklling that Cbf iiH- CTL cloud scheme 
C t ot for generalized overlap is higher than that for maximun i overfly.'- j 



When condensate heterogeneity is applied under conditions »( gwieralized o'tf&l tp (Exp. 4, 
lower right box), the effect of increased C, 0 , in the CTL cloud sc entirely cancelled out 
for CREsw through decrease in and partially -m; 


• I. ■■ out for CRS 1 1 r fl'irough i 


increase in 


F ™, . The end result is that CREsw is weaker by 2 Wnf c wiripareid t© the reference Exp. 1, 

9 

while CRElw remains stronger than in Eki> J, bui b)»n!y 1 Jw n . Note that the effect of 
inhomogeneity on CRE is stronger wit loud fraction nrv 1 . v generalized overlap (from Exp. 
3 to Exp. 4) than when i . ib**. triaxnmii:i random overlap (from Exp. 1 to Exp. 2): in the 
former case CR2$ and CREi tLemase itWreugfr/bv 6.3 Wm' 2 and 3 Wm" 2 , respectively, 
while for flic Hatter case they dd\icaie by 5 W Wm' 2 and 2.3 Wm' 2 . When the standard 
deviatic^f'uSwjfbi the beta distribucSin of condensate is halved compared to eq. (9) (box 5), 
CRE S w is reduc bout 2 WnfrAvhile CRElw is reduced by 1 Wm" 2 reflecting again the 
fact that any changes nuta fleet o^xcast fluxes instead of cloud fractions have greater impact 
on the SW compared to iS^ LW ■ 

A simultaneous change i$ both cloud fraction and condensate overlap can be achieved by 
switching from globally constant decorrelation lengths to CloudSat-based decorrelation 
lengdis (Eqs. 10-11 and Fig. 1). This process is represented by the transition from Exp. 4 to 
Exp. 8 shown by the bottom two boxes (4 and 8) of Fig. 2, CREsw strength decreases by 1 
Wm" , while CRElw decreases by 0,7 Wm" . Recall that the numbers in the boxes show 
differences with respect to the reference Exp. 1 represented by the center box (box 1), so one 
can see that transitioning from homogeneous maximum-random overlap to inhomogeneous 
clouds following a CloudSat-based generalized overlap results in 3 Wm" 2 weaker CREsw , but 
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] a slightly stronger (by 0.6 Wm 2 ) CRElw • This is possible because while cloud fraction 

2 changes (from maximum-random to generalized) have about the same effect on both the SW 

3 and LW CRE, overcast flux changes (from condensate overlap and inhomogeneity) are too 

4 weak in the LW to reverse the increased CRE of generalized overlap. 

5 The CRE response to condensate heterogeneity and generalized overlap when imposed on the 

6 cloud fields of an alternate cloud scheme can be substantially different than the one discussed 

7 above. This is shown in Fig. 3, which is the same as Fig. 2, but for the McRAS-AC cloud 

8 scheme. Cloud water inhomogeneity under conditions of maximum-random cloud fraction 

9 overlap (box 2) results in a slightly smaller weakenin&.oPw^Autf&f a slightly greater 

10 weakening of CRElw ■ The transition of homogeneous chads from ipTO, Am-random overlap 

1 1 to generalized overlap (box 3) gives a much smaller CRE ivsponfKtor McRSSiAC (4 1 Wm' 2 

12 compared to ~4 Wm' 2 for CTL). Adding inhomogenei*# W wiouds obeyu y* ig»neraJized 

13 overlap has about the same CRE effect for McRAS-AC as adding littjomogeneity to clouds 

14 following maximum-random overlap (CRE chang^-WVipfn'^xp. 3 to ifre about the same 

15 as the changes from Exp. 1 to Exp. 2); for the «TL cK-i/uwiteme -he CRE impacts diverged 

16 by 0.7 Wm' 2 ). The box corresponding!) > Trap, if mdmtcs. that when the imposed 

17 inhomogeneity is reduced by half ou'Jouds followS^jaMtalized overlap, the outcome is 

18 close to die reference effects of modified overlap and inhomogeneity 

19 largely cancel wilWIHqlbcasc umilv: CTL cloud scheme for which overlap had a 

20 much strontn;r A’RE impact than <v (Hiked inhoMPgeneity. Finally, the change from globally 

\ Jr 

21 constant declin ati on lengths to zonp lly- -dependent decorrelation lengths (Exp. 4 to Exp. 8) is 

22 notably smalle^bmhe McRAS-AC clipd fields compared to the CTL cloud scheme. 

23 This latter result is aJHiikhidcd 4 which adopts the conventions of Figs. 2 and 3, but 

24 focuses on CRE changes braujnt by changing the parameters (i.e., decorrelation lengths) of 

25 generalized overlap. The^jllt part of the figure provides global CRE impacts for the CTL 

26 cloud scheme while the right part of the figure does the same for the McRAS-AC scheme. In 

27 this figure the reference CREs come from Exp. 4 (heterogeneous clouds, generalized overlap 

28 with constant decorrelation lengths), upper left box (box 4); all other boxes contain CRE 

29 differences from these reference CREs using the sign conventions of Figs. 2 and 3. The 

30 transition from Exp. 4 to Exp. 7 (top boxes 4 and 7) captures the effect of changing the 

3 1 condensate overlap decorrelation length L r . When it is doubled from 1 to 2 km both CREsw 

32 and CRElw decrease in strength slightly. This is the result of more aligned condensate 
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distributions increasing the variability in integrated WP compared to shorter L r (more random 
overlap of layer condensate distributions producing more homogeneous WP distributions) and 
consequently yielding reduced TOAF^Jf and increased F™ . If the global decorrelation length 


of cloud fraction L a is doubled from 2 to 4 km (transition from Exp. 7 to Exp. 6, right boxes) 

the reduced C to , of the less random overlap yields further reductions of 3 Wm' 2 and 1.8 Wm' 2 
in CREsw and CRElw, respectively. Because the observed decorrelation lengths are generally 
smaller than those of Exp. 6, when they are applied in the cloud generator (transition from 
Exp. 6 to Exp. 8, bottom boxes) the CREs increase again (higher and more homogeneous 
distributions of WP) and become comparable to those of F i| 


fCTL cloud scheme, 
of the previously 


the overall impact of using CloudSat-based decorrelatk^^nglhs 

used global values of L a = 2 km and L r = 1 km (Exp. 4 to tup A . left boxes\|^wui‘,l Wm' 2 , 

slightly more for CREsw and slightly less for CRElw- Tfi£se afficiences are at'nrst glance 
rather small to justify the effort of deriving zoadly-depeaden^ ttecorpJlation lengths, 

especially since the Exp. 4 CREs are already belaxv CEWJS|FBAF and fpmore sophisticated 

M / 

treatment of overlap makes the discrepancy froifji ^serafcd c^jB^.’torse. But as will be shown 
below, the rather benign global CRE j dumges hit». local) 
substantial. /, B 


aacts that are much more 


The right part of&gp exac^Bjfk analysis as the left part, but for the McRAS- 

AC scheme j^lemented in Inii ipipact of doubling the rank correlation 

decomdatjcn&ugth (Exp. 4 to Exp^Bp about the same as for CTL, but doubling the overlap 
decorrelation 'fffltgib.does not change! RE as much for McRAS-AC. The Exp. 6 and Exp. 4 to 
Exp. 8 transitions so weaker W terms of CRE changes for McRAS-AC. When these 
results are considered m tmiunrffon with Fig. 3, the obvious conclusion is that McRAS-AC 
clouds do not cause as llg CRE changes as CTL clouds in response to the different 
prescriptions of cloud ovdrlap. We will explain why this is so in subsection 5.3. 


As a concluding thought we would like to point out that if CREsw is overestimated and CRElw 
underestimated compared to observations, as is the case for the CTL cloud scheme, it is not 
possible to bring both simultaneously closer to observations through changes in 
inhomogeneity and overlap descriptions alone. Inhomogeneity reduces CREsw and can bring 
model and observations closer, but it also reduces the already too low CRElw ■ Similarly, 
increasing CRElw via changes in overlap (i.e., increasing C to t) to match observations has the 
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1 undesired effect of making the CREsw overestimates worse. To match both components of 

2 CRE to observations, inhomogeneity and overlap changes should be accompanied by changes 

3 in other cloud properties such as cloud top height and mean condensate as well. 

4 5.2 Geographical changes in CRE 

5 In this subsection we examine whether the relatively narrow range of global CRE impact due 

6 to changes in cloud overlap specification conceals a much wider range of regional CRE 

7 changes. For the sake of brevity, we focus on only two overlap specification changes, the 

8 transition from maximum-random overlap to generalized qygrla p^yjjl t globally constant 

9 decorrelation lengths (with heterogeneous clouds), and transit: iS-Jjvtr the latter type of 

10 overlap to generalized overlap with zonally variable decstH^^fmi lengm^». parameterized 

11 per the CloudSat data analysis. In other words we examine reghjnal CRE errangr • between 

12 Exp. 2 and Exp. 4 and between Exp. 8 and Exp. 4. 

13 Fig* 5 shows maps of annually averaged CRj^sw differences belTOtr the experiments 

14 mentioned above, while Fig. 6 is a counterpari, figure Jp^MHgcp. ffhe panels in the top row 

15 correspond to Exp. 2 minus Exp. 4 diffij^fies, and «Jc pana-^i^the bottom row to Exp. 8 

16 minus Exp. 4 differences; the left pa/icl^arc for the v~fL vipud scheme and the right panels 

17 for McRAS-AC. Thc£^. .ctaudscnet^B^ds substantially greater CRE differences for the 

18 transition from u ; n-ia^Oinakio genS^fited^everlap than between two generalized 

19 overlaps, anjoti the tropics companjMb midlatpKes. Zonal CRE differences between Exp. 2 

20 and Exp* ^wik ul ~1 1 Wm‘ 2 in the »W and — 10 Win' 2 in the LW around 5°N (left panels of 

21 Fig. 7) rcflectn^Chunges in C, ot -0.13 (blue curve in the top panel of Fig. 8). The 

22 counterpart CRE diffen fotfs betwgfn Exp. 8 and Exp. 4 are ~6 Wm' 2 and —4 Wm' 2 for a C tot 

23 change of about 0.05 (re^dur. in the top panel of Fig. 8); in this case however the different 

24 vertical alignment of condensate distributions also contributes to the CRE differences, making 

25 the CREsw and CREiw changes more distinct. It is interesting that the sign of the CRE 

26 differences between Exp. 8 and Exp. 4 (changes in the details of generalized overlap) is not 

27 the same everywhere. While the CREsw ( CREiw ) difference is generally positive (negative), at 

28 midlatitudes there are negative (positive) differences with peaks at about 60 degrees latitude. 

29 The difference in behaviour from tropics to midlatitudes is solely due the parameterization of 

30 the CloudSat-based decorrelation lengths of Fig. 1. The constant decorrelation lengths are 

31 lower than those from CloudSat in the tropics and yield higher C lo , and larger WPs, ergo, 
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1 stronger CRE (expressed as positive CREsw and negative CRElw differences). In the 

2 midlatitudes on the other hand, the opposite is true, i.e., the globally constant values are above 

3 the CloudSat-based parameterized decorrelation lengths resulting in weaker CREs for Exp. 4 

4 compared to Exp. 8 (negative CREsw and positive CRElw differences). 
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The counterpart McRAS-AC CRE differences are much weaker, as can be seen in the right 
panels of Figs. 5, 6, and 7, consistent with much smaller changes in C, ot (Fig. 8) and the 
smaller global CRE differences noted earlier in Figs. 3 and 4. The zonal structure of the Exp. 
8 minus Exp. 4 CRE differences can be explained by invoking the .same arguments as earlier 
for the CTL cloud scheme, but exhibit notably smaller vaW^^e l&'H^minus Exp. 4 CRE 


for the CTL cloud scheme, but exhibit notably smaller val^^^He Is^jn^rninus Exp. 4 CRE 
differences also have the same sign as in CTL across all Rfitudes, hP^t R^l t a much weaker 
latitudinal dependence with no tropical peak as in CTL, whfL f efflg also si^Utnlially smaller. 
One interesting feature seen in the bottom panel of Fig. 'L® llmt the zonally-.; , -.raged C, oi 
difference for Exp. 2 minus Exp .4 is not only small, but genefaH«»ositive, in contrast to 
CTL. This means that either generalized overlap <bH5&B^utly results?*! ..lightly smaller total 


the zonalii 


tiall}! smaller, 
feraged C, ol 


positive, in contrast to 


cloud fractions than maximum-random overfafc. or iFrtIRte bistaiiees where C, 0 , from 
maximum-random overlap is greater exrveife trose where UH: icVerse is true. This in turn 
points to cloud vertical profiles in '’W|RAS-AC Vjfibre ih|e random part (cloudy layers 
separated by clear laym^^^^iiiri^pj^om overlap rs invoked more often than in CTL. 
Recall that withija ifee e onxSwu R ^nerah^S^/crlap, exact random cloud fraction overlap 
can only occy^i the limit of an n'Qni&lv large ® r orrelation length. 


WhyoVerlap details in the two cloud schemes affect CRE differently 


The quite distinct tire two cloud schemes when the cloud generator is 

furnished with identicairal'i t-v produce cloudy subcolumns from the same profiles of cloud 
fraction and mean condensate for radiation calculations, merits further examination. Since the 
largest impact comes from the overlap of cloud fraction, we examine here how the two 
schemes differ in terms of cloud fraction means and distributions, and the frequency of multi- 
layer cloud occurrences. 


28 First we examine the one-year cloud fraction climatology produced by the two schemes. We 

29 compare in Fig. 9 annually- and zonally-averaged cloud fraction profiles produced by CTL 

30 (top) and McRAS-AC (bottom). The differences between the two panels are striking. 

31 McRAS-AC produces in general larger cloud fractions throughout the entire extent of the 
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midlatitude and polar troposphere and the largest part of the tropical troposphere. The CTL 
cloud scheme on the other hand produces higher cloud fractions from tropical deep 
convection, and exhibits some cloud presence at the higher altitudes of the midlatitude 
atmosphere where McRAS-AC produces no clouds. The natural outcome of these average 
cloud fraction profiles is that C, 0 , is higher for the McRAS-AC cloud scheme. This is clearly 
demonstrated in the Fig. 10 zonal plot showing C !0 , from Exp. 2 (maximum-random overlap) 
and Exp. 4 (generalized overlap with L a = 2 km) which makes apparent that McRAS-AC 


produces higher zonal cloud fractions everywhere for Exp. 2 and nearly everywhere (except a 
portion of the tropics) for Exp. 4. The higher cloud fractiomWs^wh 


come with much 


greater insensitivity to the overlap specification (the distaste betwee^^ttine and red curves, 
also shown as difference in Fig. 8). Indeed, larger cloud f^a fiurfs make thsMehiils o» ; overlap 
more inconsequential since the difference between max i nwnrj . feu idom and ' my -.fegree in 
between (i.e., generalized), becomes smaller at the high end of fraction distribution. 


reduce hiultfr cloud fractions is 
tions for this quantity for 
i'nuary and July within the 
1. The seasonal differences are 
i but the differences between the two cloud 
schemes is strikir?g. McRAS-'^^^^rutes^^. t jfrore layer cloud fractions in the 0.5-0.9 
range, and, fi^ produces overcast clijud layt* which the CTL scheme never does. The 


A better way to demonstrate the tendency of Mq| 
to examine instantaneous layer cloud fractions. I 

both cloud schemes from twice-daily samples cx traced dun 

a! 

period of our runs. The four distribufi drift are shown 
not pronounced, especial® Toi^McRAva 


smaller * 1 £on>« 


t-rages of total clov 

^fe 


faction by the CTL cloud scheme in Fig. 10 appear 
therefore to be tlfc^rc^lt of consiste]^ lower than McRAS-AC occurrences of instantaneous 
layer cloud fractions' 0.5. 

Another factor that make* line details of overlap specification matter less is the number of 
cloudy layer within a gridcolumn at a particular instance. The more layers are cloudy, the 
greater the chance that they will be farther apart and therefore the greater the tendency 
towards random overlap conditions either under maximum-random overlap or generalized 
overlap. In this regard, McRAS-AC is again distinct from CTL in producing more 
occurrences of large numbers of model layers that are simultaneously cloudy (Fig. 12) within 
the gridcolumn at a particular instance. 
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1 All the above results portray a consistent picture: McRAS-AC is more cloudy than CTL under 

2 a variety of metrics and high cloud fractions are produced frequently enough to make the 

3 exact overlap specification less influential on C t0 , and CRE. 
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6 Discussion and conclusions 


global effects within climate models have been less systeratffctray qt; 
due to new capabilities in describing in within GCMs arthtrary clou# 


While earlier studies have shown that vertical cloud structure and particularly cloud fraction 
overlap can have large instantaneous effects, especially on solar flaxes (Barker et al., 1999), 

Fied. Recent progress 
lit and condensate 

overlaps that resemble more faithfully the vertical cloud ^felines observer in natule, along 
with progress on how radiation schemes can handle thes^m^c complex cloTwBPrelds, will 
help ameliorate the current state of affairs. Our study contributes to^^jneed by attempting to 
address the following question: Do the details of</^^Bto'. or lap matter ? .uji^tively to a similar 
extent when applied the exact same way on the l.taflerer 
fields produced by two distinct cloud sohn^es^The.ftlnswer] 


ctid fraction and condensate 
Igative; we find one cloud 
scheme producing cloud distribution* illicit after (uv'adaW manipulation can change the 
radiative fluxes much piqrv a n@e%(jjloy;d scheme. This means that there is no definitive 
answer on whelhecuie det*uS|pf7amid v^iOaJfcltQrcture matter much for radiation: it will 
depend on tjwrowt model and/or i’tj^wud scht^E. In contrast, introducing cloud condensate 
heteroeohcitulytound to matter moWionsistently across cloud schemes while the details of 
how the inhomoMHeo is distribution $ overlap in the vertical has only a small impact. 

The radiative conseqHawte of vertical structure and condensate heterogeneity were 
studied in this paper diagqo'ipyally, in other words, changes in radiation brought about by 
these factors did not feed J'stck to the model (a separate radiation scheme blind to our changes 
of cloud vertical correlations was running for that purpose). In that sense, our study resembles 
that of Shonk and Hogan (2010) who examined the radiative impact of different assumptions 
about condensate horizontal variability and cloud overlap by operating on cloud fields from 
re-analysis data. In that study the global effects of cloud fraction overlap (their “vertical 
shift”) on SW and LW CRE were (absolute values) ~ 4 Wm' 2 and ~2 Wm' 2 . The experiment 
transition from which these numbers were obtained are roughly equivalent to our transitions 
from Exp. 2 to Exp. 4 (see Figs. 2 and 3). In our case the change in CRE is ~3.6 Wm' 2 for 
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both the SW and LW in the CTL cloud scheme; the alternate McRAS-AC cloud scheme 
produces CRE changes slightly below 1 Wm' 2 . We conclude that studies of this type may 
eventually put an upper limit on the global impact of cloud overlap in current large scale 
models, but with a range of outcomes that may remain quite wide. Even greater variability 
range is expected to occur at smaller spatial scales. Our zonal average peak CRE impact is 
~10 Wm’ 2 , for both SW and LW CRE while that of Shonk and Hogan (2010) reaches such 
values (with much less zonal structure) only in the SW; the LW peak is about half, consistent 
with their global result. 


We did not discuss much the level of agreement of for our different 

a 

experiments with observed CRE. This was a conscious decision eement. at global 

levels at least can be achieved, through appropriate tunin'3 »u y2 and that' the best 

agreement is not necessarily achieved with the most realist!.? Miinmtions abot^lliC nature of 
the cloud fields. Nevertheless, it should be noted that if one CRE components is 

overestimated and the other underestimated, both {••mnut-he simuitarSfun’^ pushed towards 


observations by adjustments to cloud conder 
because any change that strengthens one j 
doing the same for the other components! well 
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1 Tables 


2 Table 1. Parameters for the Gaussian fits per eqs (10) and (1 1) of zonal decorrelation lengths 

3 shown in Fig. 1 . 


Fit parameters 
for eqs. (10)-(11) 

Cloud fraction 
overlap 

Condensate 

overlap 

mi 

1.43 

0.72 

m 2 

2.12 

0.79 

mj.o 

-7.00 

-8.50 | 

m 4 

-25.58 

40.40 


Table 2. List of experiments conducted with the GEOS-5 


6 schemes to assess the effects of cloud hereogeneity and overlap onx dpud radiative effect. 



inning two ooPfent cloud 


Experiment ID 


1 

MMHW ■■■L Ml 

i 2 

H 1 1 cousx lou d s 1 : ■ uux ; mum-random overlap, L r - 1 km 


Homogeneous «()|. generated overlap, L a = 2 km 

* 


5 

As 'Ivi®^. but SAth weaker cloud heterogeneity 

6 

As Exjy'4, but with L a =4 km, L r = 2 km 

7 

As Exp. 4, but with L a =2 km, L r -2 km 

8 

Heterogeneous clouds, generalized overlap from CloudSat/CALIPSO 


7 

8 
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Figure Captions 


Figure 1, (top): Cloud fraction overlap decorrelation lengths from 3° degree zonal averages of 
a(Az) for January and July 2009 (solid curves) derived from the 2B-GEOPROF-LIDAR 
CloudSat product; the dashed curves are gaussian fits according to eqs. (10) and (11). 
(bottom): As top panel, but for rank correlation decorrelation lengths calculated from 
CloudSat 2B-GEOPROF CPR reflectivities. 


Figure 2. Box chart providing diagnostic CRE changes 
distributions are changed from homogeneous to hetero- 
maximum-random to generalized for the cloud fields gei. 
scheme of GEOS-5. The changes are with respect to the r$&r 
(in Wra' 2 ) in the center box (blue for CRElw, red for C 
homogeneous clouds and maximum-random oved 
radiation package. Due to our sign convention, 
from our reference values indicate strong,!: t <JrF.', Th^|umbei 
observed values from the CERES EBAJnta set and 




clouds condensate 
changes from 
J) cloud 
'Stic CRE 
produced assuming 
ithin the RRTMG 


pbsitive CREsw changes 
'talics in the center box are 
s in parentheses are differences 


between model reference HR feign tr C HffiS observed CREs. The numbers in the left bottom 
comer of the are me expsi mustti IDs alcpi ding to Table 2. 

Figure 3. Ai.fcg. 3, but when McR' VS-kC has.yeplaced the GEOS-5 control cloud scheme. 

Figure 4. CRI- changes brought by k hanging the parameters (i.e., decorrelation lengths) of 
generalized overlap tii the CTL ylflad scheme and the McRAS-AC cloud scheme. The 
reference CREs of the >ox are from the simulation with heterogeneous clouds and 

generalized overlap with L, 2 km and L r — 1 km (Exp. 4 in Table 2). The values shown in 

the other boxes are differences from the references CREs for different experiments indicated 
by their IDs in the left bottom comer of the box according to Table 2. The left part of the 
figure is for the CTL cloud scheme while the right part for the McRAS-AC cloud scheme. 


Figure 5. Maps of annually averaged CREsw differences between the Exp. 2 and Exp. 4 (top) 
and between Exp. 8 and Exp. 4 (bottom). The left panels are for the CTL cloud scheme, while 
the right panels are for McRAS-AC. 
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Figure 6. As Fig. 5, but for CRE LW . 

Figure 7. Zonal averages of the differences shown in Figs 5 and 6. The left panels are for the 
CTL cloud scheme, the right panels are for McRAS-AC. Top panels are for CREsw , while the 
bottom panels are for CRElw- 


Figure 8. Zonally-averaged differences of C lo , (on a scale 0-100) for Exp. 2 - Exp. 4 (blue 
curves) and Exp. 8 - Exp. 4 (red curves). The top panel is for the CTL cloud scheme, while 
the bottom panel is for McRAS-AC. 


Figure 9. Annually- and zonally-averaged cloud fraction 
CTL and McRAS-AC cloud schemes. 

Figure 10. Annually- and zonally-averaged total cloud 
Exp. 2 and Exp. 4 cloud fraction overlap assumptions appl , 
schemes. 


Figure 11. Frequency distributions of twice-dai 
during January and July within the period of o 
with a separate bin for completely overca 

Figure 12. Frequency distributions 
same data as in Fig i 
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